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Abstract 

We point out that an accurate in situ determination of the earth matter density p is possible in 
neutrino factory by placing a detector at the magic baseline, L = y27c/GpN e where N e denotes 
electron number density. The accuracy of matter density determination is excellent in a region 
of relatively large #13 with fractional uncertainty 5p/p of about 0.43%, 1.3%, and < 3% at la 
CL at sin 2 20i3 = 0.1, 10~ 2 , and 3 x 10~ 3 , respectively. At smaller #13 the uncertainty depends 
upon the CP phase 5, but it remains small, 3%-7% in more than 3/4 of the entire region of 5 at 
sin 2 2#i3 = 10 -4 . The results would allow us to solve the problem of obscured CP violation due to 
the uncertainty of earth matter density in a wide range of #13 and 5. It may provide a test for the 
geophysical model of the earth, or it may serve as a method for stringent test of the MSW theory 
of neutrino propagation in matter once an accurate geophysical estimation of the matter density 
is available. 
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I. INTRODUCTION 



In the past 10 years, the atmospheric [1[, the solar [2J, and the reactor experiments [3j 
discovered neutrino oscillation and/or flavor conversion and identified neutrino masses and 
the lepton flavor mixing as the cause of the phenomena 0, HI, 0, 0] ■ The next goal of the 
neutrino oscillation experiments is to explore the yet unknown 1-3 sector of the flavor mixing 
matrix, the MNS matrix Q|, namely #13, the CP phase 5, and the neutrino mass hierarchy. 
If 6*i3 is relatively large, sin 2 20^ > 0.01, most probably conventional superbeam [9j can do 
the job. But, if it turns out to be small, sin 2 26^ 3 < 0.01, it is likely that we need a new 
technology, either the neutrino factory [Io[ or the beta beam jll|. Both options are under 



active considerations 12, 13, 14 



In measurement in neutrino factory, a long baseline distance of several thousand km to a 
detector is preferred because the event rate increases as muon beam energy gets higher and 
one wants to keep L/E not too small for oscillation signatures 15|. Then, according to the 



Mikheyev-Smirnov-Wolfenstein (MSW) theory of neutrino propagation in matter 16|, ]l7 
the earth matter effect is one of the crucial ingredients in precision measurement of the 
lepton mixing parameters. It may acts as a serious background in determination of CP 
violating phase S [T" 



191,120,121] 



In fact, it has been debated to what extent uncertainty in 
the earth matter density affects accuracies of determination of lepton mixing parameters, in 



particular, the CP phase 5. See, for example, [22], ]23|, ]24j, ]25|, [26 



One way of dealing with the problem of uncertainty of the matter density is to trust 
geophysical estimation. It appears that there is a consensus in geophysics community that 
the uncertainty in the density in the mantle region of the earth is less than 5% [281 ] . If 
the uncertainty can be regarded as robust, then the uncertainty may not be so crucial to 
hurt the accuracy of measurement of the mixing parameters, though it still produces sizable 
uncertainties. Notice, however, that the authors of [24J take much more conservative attitude 



in estimating the uncertainty of the matter density. See e.g., [27| for a recent estimation of 
impact of the matter density uncertainty on accuracy in mixing parameter determination. 

However, it would be much nicer if neutrino factory experiments can measure in situ the 
average matter density in the mantle region of the earth. It is not completely satisfactory 
that the most advanced apparatus for precision measurement of the lepton mixing param- 
eters has to rely on the parameter that cannot be directly measured. In fact, the same 
attitude is shared by the authors of the pioneering work of neutrino factory analysis (l5| . 
in which the accuracy of matter density determination by neutrino factory is estimated. 
Describing an improved way for the in situ measurement of the earth matter density in 
neutrino factory is nothing but the goal of this paper. The recent progress in understanding 
of how to lower the analysis energy threshold 2^] will be the key to the remarkable accuracy 
we will uncover. 

Moreover, measured matter density in deep interior of the earth by neutrinos would offer 
an independent test for the theory of earth structure and its formation. Note that the 
mantle region contains 70% of the earth mass. If it is significantly different from the value 
predicted by geophysical models (after taking account of the uncertainties in both prediction 
and measurement), we have to think about at least one of the following possibilities: (1) 
A drastic revision of geophysical model of the earth must be attempted. (2) The MSW 
theory of neutrino propagation is in error. We hope that this discussion is illuminative of 
the importance of the neutrino measurement of the earth matter density. 

The chondrite earth model with constraint of seismic data results in earth model with 
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several layers and the matter density in each layers can be predicted with reasonable precision 
with constraint of accurately measured earth mass. In particular, the density in the mantle 
region seems to be obtained with the least uncertainties [30]. If the geophysical estimation 
can be regarded as robust neutrino factory measurement of matter density we propose in this 
paper severely tests the theory of neutrino propagation in matter. Notice that the matter 
density at the solar core is currently measured by using neutrinos only up to a factor of ~2 
uncertainty |3l| . Neutrino tomography of the earth by using accelerator and other neutrino 



sources has been discussed in many literatures [32], [33|, [3J, |35[ . 

In Sec. HIl we give three conditions for accurate measurement of the earth matter density 
based on general considerations. In Sec. IHIt we discuss what is the baseline that is most 
sensitive to matter density change. In Sec. HVt we describe a concrete way of determining 
the earth matter density by neutrino factory. In Sec. |V] we define the statistical method 
for analysis. In Sec. |VI| we describe the results of our analysis. In Sec. IVII[ we address 
the problem of 5-dependence of the error of matter density determination. In Sec. IVIII[ an 
iterative procedure for combined analysis of data at the intermediate and the far detectors 
is described. In Sec. IIX[ we give concluding remarks. 



II. GENERAL CONSIDERATIONS 

Experiments which aim at measuring the earth matter density accurately by using neu- 
trino oscillation must satisfy a few crucial requirements. Firstly, the matter effect must be 
sizable, at least comparable with the vacuum oscillation effect. The relative importance 
between the vacuum and the matter effects in long-baseline (LBL) neutrino oscillation ex- 
periments may be quantified by comparing the following two dimensionless quantities: 

\Am 2 31 \L _ f\Am 2 31 \\f L \ ( E V 1 
Asi = ~^E—- l21 \W^) UoOOkmV (jGevJ U) 

aL = \G F N e L = 0.27 (—^s) (tJ—) (2) 



y/2 ' V2.8g/cmV ViOOOkm 

where a is related to neutrino's index of refraction with Gp being the Fermi constant and 
iV e the electron number density in the earth. N e is related to the matter density p as 
N e = Y £ NaP with Avogadro's number Na and the electron fraction Y e . We take Y e = 0.5. 
In most cases, the experimental setup is such that A 31 ~ 7r/2 to maximize the appearance 
signal. In view of (j2j), the matter effect plays no significant role in settings of baseline less 
than ~ 1000 km. 

Sizable matter effect is certainly the necessary but not the sufficient condition for ex- 
periments to measure the matter density. What is important is the sensitivity to artificial 
matter density change Sp. To quantify the sensitivity one has to have a variable experimen- 
tal parameters through which one can control relative importance of the matter effect to 
that of the vacuum oscillation. Since aL/A 3 i is proportional to E the natural candidate for 
relevant experimental parameters to vary would be the neutrino energy. 

It is not the end of the list of requirements. A measurement at a given set up determines 
a combination of the mixing parameters, typically #13, S, and the matter effect coefficient 
a. Unless the former two parameters are separately measured in experiments in vacuum, 1 

1 MEMPHYS [36| would be an ideal apparatus for this purpose because even the T2K experiment [37| is 
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the coefficient a or the matter density alone cannot be cleanly determined by any single 
measurement. Generally speaking, the condition is hard to be met because most of the 
projects to determine unknowns in the lepton flavor mixing are designed to be sensitive not 
only to the matter effect, which is required to determine the mass hierarchy, but also to the 
CP phase 8. 

In summary, we list here the three requirements (A) — (C); (A) sensitivity to the matter 
effect. (B) existence of dp sensitive tunable parameters. (C) mixing parameter independent 
measurement of the matter density. We will discuss in the following sections below how 
these requirements can be satisfied. 



III. WHICH BASELINE?; PRELIMINARY CONSIDERATION WITH ENERGY 
SCAN 

In this paper, we focus on neutrino factory measurement of the earth matter density. 
Since most part of the neutrino trajectory in neutrino factory is in the mantle region of 
the earth, what we mean by the matter density in this paper is that of the mantle region. 
Precisely speaking, the matter density we try to determine is an averaged value along the 
neutrino trajectory. 

The first question we must address is; What is the appropriate baseline distance for 
the measurement? We try to find an answer to this question by choosing muon beam 
energy as the tunable parameter, the energy scan. Namely, we consider the neutrino factory 
measurement of numbers of appearance events v e — > at two different muon beam energies. 
It will give us the normalized event number difference AN/N. The question we address here 
is what is the baseline L that gives rise to the strongest response in AN/N to matter density 
change. 

For simplicity, we employ in this section the asymptotic expansion of the oscillation 
probability assuming A 31 C l. 2 The standard formula for v e — » appearance probability 
[15( (for its full expression, see f[T2l in Sec. IVIII) takes the form at high energies, A 31 = 



^ ] « I . as 



A B 

P(^^) = e5 + W (3) 

where E is the neutrino energy and L is the baseline distance. The coefficients A and B are 
given to leading order in s 13 by 

(sin clL "\ ^ f \ /\ 777, ^ I j^j \ ^ 
~^L~) [ ~AJ [ 4s 23 s ?3 ±4sin2^ 12 c 2 3S 2 3es 13 cos5 + ^3sin 2 20 12 e 2 ] 



, sin aL\ 2 ( I Ami \L\ 3 . . 



2s 23 s i3 + sin 26 l i 2 c 23 s 23 esi 3 cos <5 ( I ± 



tan S 
g(aL) 



(4) 



where e = | A m jP |, and a = GpN e /y/2 as before. The function g is defined by g(x) = 
1/x — cot x. The positive and the negative sign in (j3J) are for the normal and the inverted 



contaminated by the matter effect which leads to some sensitivities to the mass hierarchy in a limited 



region of 5, as anticipated 2l|, |38( and proved by the recent analyses [39j |. 
^ In a standard setting, A 31 is small; A 31 = 0.24 ( ^fj^ ) (aoofflEBi) (wfev) 
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mass hierarchies, respectively. In (j4]) and throughout this paper, we use the constant matter 
density approximation. 

It is well known that the integrated electron neutrino flux F ve from muon decay has 
energy distribution as F ue (E) = F y 2 (l — y) (in units of number of z/'s per unit area) where 
y = and F = ^ jjs§$ with and L being the number of useful muon decay and the 
baseline distance, respectively. Then, the number of wrong sign muon events is given, using 
the number of target atom Nt, and by approximating the energy dependence of the charged 
current (CC) cross section as linear, <j C c — &oE, as 

N(E^) = N T / dEF e (E)a C c(E)P(is e -> ^, E) = -—^a N T (E^A + SB) (5) 

/ W , IX Li TTl,, 

where we have assumed that the muon's threshold energy E t h can be made low enough to 
allow the approximation 1 S> ^ — 0. 

Let us consider measurement at two adjacent muon beam energies E^ and E^ + AE^ with 
number of events N(Efj) and N(E^ + AE^), and define their difference as AiV = N(E^ + 
AE^) - N(E^) ~ ^§^AE^. We define the double ratio ^/^, which is independent of 
the uncertainties in the CC cross section, the neutrino flex, the baseline, and the number of 
target atoms. We obtain 



AN_ 1 



AE„ , , o_B_ 



(6) 



Now, the crucial question we have to ask is: What is the baseline distance for which 
variation of the matter density, or a, induces maximal changes in the double ratio? From 
(EJ), it is the place where B/A changes significantly when aL is varied. The answer is aL = n 
at which g{aL) (which is proportional to B/A) diverges. This distance is nothing but the 
one called the "magic baseline" in the literature [4(|. It is interesting to see that the magic 
baseline appears in our treatment as a point where the sensitivity to matter density variation 
is maximal. The length has been known in the theory of neutrino propagation in matter as 



"refraction length" [16j |. For a recent discussion on the meaning of the magic baseline, see 

Despite that the asymptotic expansion is not valid at around the distance where g(aL) 
diverges, one can show by using the full expression of the oscillation probability that the 
distance comparable to the magic baseline is indeed the most sensitive place for the double 
ratio to change in matter density. (See Fig. 5.1 of [43] .) Thus, the asymptotic expansion 
is a correct indicator for the right baseline distance for measurement of the matter density. 
(It is reminiscent of the feature that the one-loop QCD coupling constant diverges at the 
hadronic scale A.) 



IV. MEASURING EARTH MATTER DENSITY AT THE MAGIC BASELINE; 
ENERGY BINNING 

We now switch to a different strategy of using neutrino energy as the tunable parameter, 
though we will make a brief comment on the method of energy scan at the end of Sec. IVIl 
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We consider measurement of the earth matter density at the magic baseline 



V2- 



71 



G F N P 



7480 



P 



4.2g/cnT 



km. 



(7) 



As was emphasized in the original article [40j, one of the most characteristic feature of 
the magic baseline is that the oscillation probability P[y e — > v^) is independent of the CP 
violating phase S [43[ . Then, one can measure #13 independently of 5, and this property 
has been utilized to resolve the parameter degeneracy by combining with detector at L — 
3000 — 4000 km in neutrino factory measurement of #13 and 5 We will show below 

that the property is also the key to our method for accurate measurement of the earth 
matter density, fulfilling the requirement (C) at least partly. We, however, encounter below 
the problem of unexpected ^-dependence in a limited region of # 13 and 5. See Sec. IVHI 

Now, the key question is: What is the most efficient way to measure the matter density 
accurately at the magic baseline? To formulate the right strategy for this purpose we 
analyze the structure of the appearance probability P[y e —>■ v^). At around the magic 
baseline the Cervera et al. formula 15] has a very simple form with vanishingly small solar 
and interference terms (the first term in ( fl2l) in Sec. IVHI) . 



S23 sin 2 29 



13 



A31 sin(aL =F A31) 
(aL =F A31) 



(8) 



with aL = 7r, where the ± sign in ([8]) corresponds to the neutrino and the anti- neutrino 
channels. We note that when the matter density is perturbed, p — > p + Sp, or aL = 71 + e 
at the magic baseline, the response of the function in the square bracket in (jSj) is given as 
follows: 



sin(aL =f A 31 ) _ sin(7r =F A31) 
(aL =f A 3 i) ~ (71-T-A31) 



[1 - eg(7r T A 3 i)] 



(9) 



where g(x) = 1/x — cotx as before. 3 Equation ([9]) indicates that the response to density 
change depends upon the neutrino energy, the channel [y or u), and the neutrino mass 
ordering, the normal (A31 > 0), or the inverted (A31 < 0) hierarchies. 

For definiteness, we first discuss the case of the normal hierarchy (A31 > 0). In the 
neutrino channel, the equation g(n — A 31 ) = has a solution at A 31 — 7r = 0. Namely, 
the critical energy is given by E c ~ 7.6 GeV for L = 7500 km. Since g(x) is an odd 
function of x, the probability decreases (increases) at neutrino energies E > E c (E < E c ) 
as the density increases in the neutrino channel. In the anti-neutrino channel, the equation 
g(n + A31) = has a solution at A31 + tt = 4.493, and g(x) > for x > 4.493 and vice 
versa. The corresponding critical energy is E c ~ 17.6 GeV. The antineutrino probability 
increases (decreases) at neutrino energies E > E c (E < E c ) as the density increases, a 
behavior opposite to the neutrino channel. In the inverted hierarchy, the situation of the 
neutrino and the anti-neutrino channels completely reverses, as one can easily confirm. 



It may be instructive to remark that g(x) is a monotonically increasing function of a; in a range 2mr < 
x < (2n + 1)tt where n is an integer. It has zero at x = and at x = 4.493, and diverges to ±00 at 
x = nir =F e. 
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FIG. 1: The energy distribution of event number (per GeV) is plotted with three values of the 
matter density, 4.2 g/cm 3 (shown in blue dotted curve), 4.3 g/cm 3 (black solid), and 4.4 g/cm 3 
(red dash-dotted). The left and the right two panels in Fig. Q] are for the cases of the normal 
and the inverted mass hierarchies, respectively. The mixing parameters are taken as 5 = and 
sin 2 26»i3 = 0.01. 

In Fig. [TJ the event number distributions are plotted as a function of neutrino energy 
with three values of the matter densities. 4.2, 4.3, and 4.4 g/cm 3 . The mixing parameters 
are taken as 5 = and sin 2 2^13 = 0.01. The left and the right two panels of Figs. [Hare 
for the cases of the normal and the inverted mass hierarchies, respectively. Notice that the 
oscillation probability is the only piece that is sensitive to matter density variation, and hence 
our above estimation of the critical energy should apply to the event number distribution as 
well. Figure [1] confirms the qualitative behavior expected by the above analysis of response 
of the oscillation probability to the matter density perturbation. The difference between the 
critical energies in the neutrino and the anti-neutrino channel, as well as the characteristic 
feature is obvious; Smaller (larger) number of events in the neutrino (anti-neutrino) channel 
at energies above E c for higher matter density in the normal hierarchy case (and a completely 
reversed feature in the inverted hierarchy case). Therefore, it appears that one can expect 
the largest sensitivity to the matter density by using the two energy bins below and above 
E c . The opposite response of the neutrino and the antineutrino channels to the matter 
density variation suggests that higher sensitivity can be achieved by using both channels, 
which will be verified in Sec. IVII 

We emphasize that lowering the analysis threshold, the possibility under active inves- 
tigation 29[, is the key to the feasibility of our method of energy binning, because the 



information below E < 10 GeV is quite essential for neutrino (anti-neutrino) channel for 
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the normal (inverted) hierarchy. We used the formula (jSJ) but with the exact expression of 
the oscillation probability [44] to compute the number of events with suitable modification 
of the range of integration over neutrino energy. In our analysis in this paper we use only 
the signal events by the charged current (CC) reaction. We take the cross section of CC 
reactions (including deep inelastic reactions) occ — 0.67 x 10~ 38 (£'/GeV) cm 2 for neutrinos 
and ace — 0-34 x 10 _38 (-E'/GeV) cm 2 for anti-neutrinos [15) in our analysis. Throughout 
this paper, apart from the comment at the end of Sec. IVI Al the true values of the neutrino 
mixing parameters are assumed as follows: jAm^l = 2.5 x 10~ 3 eV 2 , Am 2 x = 7.9 x lCT 5 eV 2 , 
sin 2 29 12 = 0.86, and sin 2 2 9 23 = 1. 

In principle the disappearance channel — > z/ M can be added to the analysis within the 
setting of magnetized iron detector. We do not go into this task in this paper; Though the 
matter effect is known to be sizable at these distances [45j the response to matter density 
change does not appear to be favorable. 



V. ANALYSIS METHOD 
A. Experimental setting 

In this paper, we consider a setting of neutrino factory with muon energy E., = 50 GeV 
and assume total 3 x 10 21 useful muon decays per each polarity, as quoted in j46|. It may 
require 5-10 years operation of muon storage ring depending upon its luminosity. 4 Our basic 
attitude is that operation of the neutrino factory will be coordinated by the requirement of 
optimization for accurate measurement of CP phase 5 and #13. Therefore, measurement of 
the matter density should be carried out as a "byproduct" , within a given schedule optimized 
for the above purpose. 

A standard setting of the neutrino factory includes two magnetized iron detectors, one at 
baseline L = 3000 — 4000 km and the other at the magic baseline L = 7500 km, which will be 
denoted as the intermediate and the far detectors, respectively. Note that the distance to the 
far detector does not exactly coincide with the magic baseline as given in (j7j) if the matter 
density is 4.3g/cm 2 , for example. But, the discrepancy of this level should be admitted 
because it is unlikely that we can find a site for the far detector at the distance equal to the 
magic baseline in a mathematical precision. Note also that for L = 7500 km the maximum 
depth of the neutrino trajectory is about 1220 km in the lower mantle region of the earth. 
Hence, neutrino passes through the lower mantle region in more than half a fraction of its 
trajectory. Note also that for L = 7500 km neutrino passes through the lower mantle region 
in more than half a fraction of its trajectory and its maximum depth is about 1220 km. 

In this paper, we analyze only the data taken by the far detector, leaving a full analysis 
of the data at the intermediate and the far detectors for future work. As a way of effec- 
tively implementing the information obtained by the intermediate detector, we impose the 
constraint on 5 in doing fit in the analysis of far detector data. The fiducial mass of the far 
detector is assumed to be 40 kton. We also assume that the efficiency of event reconstruc- 
tion is 80% and is independent of energy, which appears to be a good approximation to the 



The authors of 25 



27| assume more modest value of total 2 x 10 21 useful muon decays per each polarity. 
On the other hand, if the "reference neutrino factory" setting [47J is realized one may be able to accumulate 
5 x 10 21 muon decays by 5 years operation per each polarity. 
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efficiency in most of the relevant energy region shown in 2{|. 5 See also 48] for the plot. 

We now argue that analyzing only the data of the far detector is a sensible first approx- 
imation. To a good approximation there is a separation between the intermediate and the 
far detectors about their functions; The far detector has no sensitivity to the CP violating 
phase 5, and the sensitivity to it solely comes from the intermediate detector. The matter 
effect which is crucial to determine the earth matter density is mainly felt by the far detec- 
tor. Though 0i3 is determined jointly by the intermediate and the far detectors, the latter is 
crucial to resolve degeneracy to precisely measure #13 and 5 [40]. In fact, we will show that 
the maximally achievable accuracy of determination of the matter density only with the far 
detector is already quite remarkable, ~1% (~2%) level for sin 2 2#i3 = 0.01 (0.001) at 1 a 
CL, as will be shown in Sec. IVII 



B. Statistical procedure for quantitative analysis 

In our analysis, we use two energy bins, 5 GeV < E < 10 GeV and 10 GeV < E < 50 GeV 
for neutrino, and 5 GeV < E < 20 GeV and 20 GeV < E < 50 GeV for anti-neutrino 
channels for the reasons discussed in Sec. OV] Notice that the boundary between the low- 
and the high-energy bins are somewhat different between the neutrino and the anti-neutrino 
channels, as can be seen to be appropriate in Fig. [TJ While we do not explicitly deal with 
the issue of finite resolution in the reconstructed neutrino energy it will be partially taken 
care of by the systematic error discussed below. 

The definition of Ax 2 for our analysis is given by 



Ax 2 = min > 



2 2 1 ^,2 



V j W - (1 + Qt + Q a + a)K? P ) of I *^ 
^ 1 NT + ^ b {NT) 2 *M -Is 



a 2 



+ -5-,(10) 



a PB 



where N°^ s is the number of events in z-th bin computed with the values of input parameters, 
and N^ p is the one computed with certain trial set of parameters. The form of Ax 2 with 
the notation of errors is analogous to that used in [49||. The generalization to measurement 
at multiple energy bins, if necessary, is straightforward. 

The nature of the errors <Jpb, crpb, &pB, arid a p b and the examples of them are as follows. 
For simplicity of notation, we call the following errors as the category A; detection efficiency, 
energy resolution, and efficiency in muon charge identification. We assume that the dominant 
part of the category A errors are correlated between v and v channels. This assumption can 
be relaxed, if necessary. 

• a PB [y — v and bin-by-bin correlated error): Uncertainties in detector volume, energy- 
independent component of the category A errors. 

• <^Pb [y — v correlated but bin-by-bin uncorrelated error): Energy- dependent component 
of the category A errors. 



5 The efficiency presented in these references is a preliminary one and is being under further improvement. 
The currently estimated efficiency may be lower than 80% at energies at around 5 GeV. However, our 
optimistic assumption may not hurt the sensitivity estimate significantly because the event numbers which 
comes from this region is rather small. 



9 



• crpB [y — v uncorrelated but bin- by-bin correlated error): Uncertainties in beam energy 
estimate, neutrino flux estimate, energy-independent part of reaction cross section 
error. 

• <jpb (uncorrelated error): Energy- and channel-dependent fluctuation of the category 
A errors, etc., which should be very small. 

We take, without any solid informations for these errors at this moment, opb = &Pb = 
a P B = 2% and a p b = 1% in our analysis. The expected small values of the errors are 
partly based on the fact that the uncertainties in the muon energy and the luminosity are 



negligibly small for our present purpose [12j, |50|. We also examine the stability of the results 



by relaxing these errors by a factor of 2. We will see that measurement is still dominated 
by the statistical error and a much smaller uncorrelated error, a pb = 0.1%, would lead to a 
minor improvement. 



C. Treatment of CP phase 6 

Despite that ^-dependence is expected to be absent at the magic baseline the proper 
procedure of the analysis has to include varying over 5 during the fit. The best way to do 
this is to carry out a combined analysis of data taken by both the intermediate and the 
far detectors. Leaving it to a future work, we give in this paper a simplified treatment by 
mimicking data at the intermediate detector by a constraint imposed on S by adding the 
term 

A^=(^|^) 2 . (11) 

in Ax 2 in (fTOl) . We have chosen as = 0.35 (20 degree) as a representative value based on 
the estimation of the width of \ 2 parabola at around the local minimum as given in 51 



We note that "by definition" of the magic baseline the procedure of varying over 5 should 
give no significant effect. In fact, we have confirmed that the property holds in a wide region 
of #i3 and 5, apart from the one with a subtle feature to be addressed in Sec. IVHI We have 
also confirmed that using a factor of 2 larger value of as does not produce visible changes 
in the results in the entire region. 



VI. ANALYSIS RESULTS 
A. Analysis with fixed 5 

We first discuss the results with fixed 8. It is to have better understanding of the structure 
of analysis of the far detector and the roles of the systematic errors without worrying about 
influence of 5. 

In Fig. [2], presented is the region in sin 2 29% 3 — p space allowed by measurement of and 
i>fj, appearance events with our reference setting of neutrino factory. The confidence level 
(CL) contours are defined with 2 degrees of freedom (DOF) in Fig. [2j The input matter 
density is taken as p = 4.3 g/cm 3 . The normal mass hierarchy is assumed and 5 is taken 
to be zero. The true values of sin 2 26*i 3 is assumed to be 0.01 and 0.001 in the left and 
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FIG. 2: Regions in sin 2 2#i3 — p space allowed by operation of neutrino factory with 3 x 10 21 useful 
muon decays per each polarity watched by 40 kton magnetized iron detector at L = 7500 km. The 
red solid, the green dash-dotted, and the blue dotted lines are for contours at 1, 2, and 3 a CL, 
respectively, defined with 2 DOF. The left (right) two panels are for sin 2 2#i3 = 0.01 (0.001). The 
upper two panels are for the systematic errors aps = opb = &pB = 2% and a p b = 1% as described 
in the text, and the lower two panels are with a factor of 2 larger (apply to all) errors. The normal 
mass hierarchy is assumed and CP phase is taken to be 5 = 0. 

right two panels, respectively. The upper and the lower two panels are for the systematic 
errors given in Sec. IV Bl and the twice larger values than those, respectively. Notice that 
the region is compact not only in p direction but also in sin 2 2# 13 direction despite the fact 
that we have analyzed the data at the far detector only. (If the binning of the data were 
not done the analysis would yield a prolonged contour in sin 2 29i 3 — p space, as one expects 
from ([8]).) It is due to a high sensitivity to # 13 at the magic baseline which stems from that 
the probability is free from 5. It is impressive to see this behavior because the statistics is 
quite small due to the long baseline, as can be seen in Fig. [U 

As indicated in the lower two panels in Fig. [21 a factor of 2 enlarged systematic errors gives 
a relatively minor effect on Sp/p, making the original uncertainty 1.3% (2.8%) to a modestly 
larger value 1.4% both at 1 a CL at sin 2 2#i 3 = 0.01 (0.001). On the other hand, 

the enlarged systematic errors gives a larger effect on #13, worsening the original uncertainty 
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of sin 2 20i3 from (±1%) to (+} 2 %) to (+ 2 ^%)] at 1 a CL at sin 2 20 i3 = 0.01 [0.001]. 

The fact that a factor of 2 relaxed systematic errors lead to only modest increase of Sp/p 
must be emphasized. The tendency is more prominent for sin 2 26 ri = 0.001. It is because 
the detector at the magic baseline is highly sensitive to density change and measurement is 
still dominated by the statistical error. 

Different roles played by various errors which lead to the enhanced uncertainties should 
be mentioned. The factor 2 larger overall normalization error opb primarily affects the 
uncertainty in sin 2 26 l 13 as expected, and gives little effect to the uncertainty of matter 
density 8 p. The increasing error of Sp/p is due to accumulation of small and comparable 
contributions by the other three kind of the systematic errors, apt,, o v b and a p \j. The 
enhanced error of sin 2 2#i 3 , on the other hand, is contributed mostly by Opb and a pE apart 
from minor contribution from <jpb- 

Here are some comments on the effects of including the energy resolution and possible 
values of Arn^j different from the standard one. The numbers presented below are for the 
case of normal mass hierarchy but we have checked that the situation is similar in the case 
of inverted hierarchy. The CP phase S is taken to be 5 = 0. 

(1) Because we use only two wide ranged energy bins, it is likely that including the energy 
resolution into our analysis gives only a minor effect to the sensitivity to Sp/p. To confirm 
this point explicitly we have included the energy resolution into the analysis by assuming 
the gaussian form with width <je = 0.15E GeV (as taken in [271] )■ At sin 2 28 13 = 0.01, 
Sp/p becomes worse from 0.81% (without erg) to 0.86% (with a E ). At sin 2 26*i 3 = 0.001, the 
corresponding numbers are 1.8% (without <te) an d 1.9% (with erg). Since the effect of <je is 
so minor we ignore the energy resolution in our subsequent analysis. 

(2) If the value of Am|i is larger than Arn^j = 2.5 x 10 -3 eV 2 which is assumed in our 
analysis, there is no problem; Sp/ p is comparable or smaller than the one obtained above. 
However, If the value of Am 31 is smaller, the sensitivity to Sp/ p becomes less. If we assume 
Amgj = 1.5 x 10 -3 eV 2 , the lower limit allowed by the Super-Kamiokande data (the last 
reference in |lj), we have to lower the critical energy E c in defining the energy bins. We 
have used E c = 8 GeV and E c — 10 GeV for the neutrino and the antineutrino channels, 
respectively. The obtained sensitivity to Sp/p becomes worse but only moderately from 
the ones quoted above to 1% and 2% at sin 2 2#i 3 = 0.01 and 0.001, respectively. We note, 
however, that the latest value of Am^ from the MINOS experiment is higher, Am^ = 
2.74±q;25 eV 2 (la CL), and that the accuracy of determination must be greatly improved 
by the T2K and the NOz/A experiments. 

B. Analysis with varied 5 

We now turn to the analysis with varied S; During the fit we vary S with addition of 
Axj i n (JHJ) . In Fig. [31 we present the fractional error Sp/p as a function of sin 2 2# 13 at 
1, 2, and 3 a CL. Hereafter, the uncertainty Sp is defined with 1 DOF by marginalizing 
sin 2 2#i 3 . The upper and the lower two panels in Fig. [3] are for the normal mass hierarchy 
with S = and for the inverted mass hierarchy with S = 47r/3, respectively. As we will fully 
discuss in Sec. IVHI the accuracy is the best at around S = and S = it in the normal and 
the inverted hierarchies, respectively. The case of S = 47r/3 is shown in Fig. [3] to indicate 
that the 5-dependence (which exists only in small #i 3 region) is mild and to avoid to confine 
ourselves into the CP conserving values of S. It should be mentioned here that though we 
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sin 20 13 



FIG. 3: The fractional errors in the matter density determination dp/p at 1, 2, and 3 a CL defined 
with 1 DOF by marginalizing Bis are plotted as a function of sin 2 26±3 by the red solid, the green 
dash-dotted, and the blue dotted lines, respectively. The upper panel is for the normal mass 
hierarchy with 5 = and the lower panel for the inverted mass hierarchy with 5 = 47r/3. 
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have varied S for proper analysis procedure it does not produce, apart from region of very 
small #13, any sizable changes in the results, in accord with the natural expectation. 

Notice the remarkable accuracy of determination of the matter density p at the magic 
baseline represented in Fig. [3J The uncertainty is only about 1% at sin 2 29% 3 = 0.01 at la CL 
for both the normal and the inverted mass hierarchies. At sin 2 2# 13 = 0.001, the uncertainty 
increases to about 2% (2.5%) at la CL for the case of the normal (inverted) mass hierarchy. 
At sin 2 26*13 — 0.0001, however, Sp/p becomes worse to about 3% (4%) at la CL for the case 
of the normal (inverted) mass hierarchy. 

Unfortunately, however, we have to mention that a subtlety exists on ^-dependence of 
the sensitivity, which is quite unexpected from the usually advertised nature of the magic 
baseline. It will be fully discussed in the next section. 

Before entering into the problem, let us note the followings: 

(1) We have explicitly confirmed by plotting Sp/p as a function of L that distances compa- 
rable to the magic baseline are the right distances for accurate determination of the earth 
matter density in the energy binning method. In fact, the fractional error \Sp/p\ has broad 
minima in region of distance L = 7000 — 9000 km if S is fixed, and L = 7500 — 9000 km 



if S is varied. See Fig. 5.5 in 42|. (The former exercise is to avoid extra complication due 
to S dependence which exists in the oscillation probabilities in general except for the magic 
baseline.) 

(2) We have also investigated the possibility that the accuracy of determination of the earth 
matter density can be improved by doing energy scan as well as neutrino energy binning, 
but without success. 



VII. 5 DEPENDENCE OF UNCERTAINTY IN MATTER DENSITY DETERMI- 
NATION 

Despite the remarkable accuracy presented in Fig. [3] it is not the end of our work. We 
must address the issue of a rather strong (^-dependence in the fractional error Sp/p at small 
#13, which significantly alters the feature of the results presented above, as the "best cases" 
with the particular values of S chosen. 

In Fig. HI presented is the fractional uncertainty Sp/p at 1 a CL defined with 1 DOF as a 
function of 5 for four different values of 6*13, sin 2 29 13 = 0.1 (red solid line), sin 2 26' 13 = 0.01 
(green dashed line), sin 2 26*i 3 = 0.003 (blue short-dashed line), sin 2 26*i 3 = 0.001 (magenta 
dotted line), sin 2 29 13 = 0.0001 (light-blue dash-dotted line). There is a significant differences 
between the cases of small and large #13 separated by a critical value, sin 2 29 13 ~ a fewx 10~ 3 . 
At larger 9\ 3 than the critical value, the sensitivity to matter density is excellent; The 
fractional error Sp/p is about 0.43%, 1.3%, and < 3% at la CL at sin 2 29 l3 = 0.1, 10" 2 , 
and 3 x 10~ 3 , respectively, for all values of S. It is also notable that (^-dependence of Sp/p is 
very mild in the large 9i 3 region, which confirms the naive expectation due to independence 
on S at the magic baseline. 

At small #13, however, there exists a significant 5-dependence in Sp/p. Despite that the 
Sp/ p remains rather small in most of the region of S, about 4%— 6% even at sin 2 2# 13 = 10~ 4 , 
there are some "spike" structures in Sp/p at around S ~ tt for the normal hierarchy and 
S ~ for the inverted hierarchy around which Sp/p blows up to 15%— 20%. They arise when 
a separated "island" that appear in la CL equi-% 2 contour merges with the "mainland" 
la CL allowed region. In our analysis we do not make a sophisticated treatment to take 
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into account separated islands in estimating the errors. Therefore, 5p/ p in regions near the 
spike is not quantitatively reliable. On the other hand, the intricate structure of the equi- 
X 2 contour suggest that it is unstable to inclusion of additional informations. Namely, a 
significant improvement of the sensitivity can be expected in small #13 region by a combined 
analysis of the intermediate and the far detectors. 

Apart from the small structures there seems to be a relationship between the values of 5 
with worst sensitivity in the normal and the inverted mass hierarchies as <5 (worst ) norma ; m 
5 (worst) inverted +7T- The strong ^-dependence of the sensitivity including this feature is quite 
an unexpected subtlety because we usually think of the magic baseline as the distance where 
the CP phase 5 plays minor role. Therefore, the feature should be understood better and it 
is the reason why we devote the rest of the present section to analyzing the subtlety. 



The standard form of oscillation probabilities in 
by 



v e and v a — > v e channels are given 
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where ± indicates the sig n of Am| x and A 21 = Am^L/AE. In the last line in fll3p we have 
used the relationship [19| between the coefficients X± and Y± in neutrino and antineutrino 
channels, X± = X^ and Y± = —Y T . 

We perturb the matter density around the one corresponding to the magic baseline as 
aL = 7r + e. We expand the neutrino and the antineutrino oscillation probabilities by e and 
obtain and obtain, to first order in e, 
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Suppose that values of the parameters are such that the quantity in the square bracket in ([Tj 
or (fT5|) cancel out. Then, the measurement loses the sensitivity to density change, leading 
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FIG. 4: Presented are the fractional errors Sp/p at 1 a CL with 1 DOF as a function of 5 for 
five different values of #13, sin 2 29\3 = 0.1 (red solid line), sin 2 2^13 = 0.01 (green dashed line), 
sin 2 26*13 = 0.003 (blue short-dashed line), sin 2 26>i 3 = 0.001 (magenta dotted line), sin 2 26>i 3 = 
0.0001 (light-blue dash-dotted line). The upper and the lower panels are for the normal and the 
inverted mass hierarchies, respectively. 
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to an enlarged error in 5p/ p. The phenomenon can take place only if the two terms have the 
same order of magnitudes, which occurs at small #13, sin 2#i3 m A21/A31 = \Am\ l / Am 2 ^ ~ 
1/30, which explains why the (^-dependent loss of the sensitivity in Fig. H] occurs only in 
region of sin 2 26 13 comparable to ~ 0.001. 

We note that the event statistics is dominated by the neutrino channel in the normal mass 
hierarchy, and by the antineutrino channel in the inverted one, one can observe in Fig. [TJ 
Let us focus on these dominant channels. Then, we observe a notable regularity; The values 
of S which makes the square brackets in (fTTj) (upper sign) and (fT8j) (lower sign) vanish are 
related by cos(5 l/;norma i - A 3i ) + cos(5 P;inv erted - A 31 ) = 0, or 8 u;noimai = 5p- invevted + ir. It 
explains the salient feature of Fig. Hf, the relationship between the 5's corresponding to the 
worst sensitivities to 5p/p in the normal and the inverted mass hierarchies. One can also 
verify that with A 31 ~ 1 the equations of vanishing the order e terms in (fTTj) and (fTSl) are 
roughly consistent with 8 u . >msrvaai ~ 7r and S U]invCTtcd « 0. 



VIII. GLOBAL FIT VS. ITERATIVE METHOD 

In this paper, we have restricted ourselves to the simplified analysis in which information 
from the intermediate detector is modeled by the effective \ 2 as in ffTTl) . Clearly one has to 
engage in the next step a global analysis of all the data set taken by the intermediate and 
the far detectors which aims at measuring the earth matter density in situ as a part of the 
program of determining all the relevant lepton flavor mixing parameters. 

What is the right way to carry it out? A straightforward answer to this question may be 
to marginalize (minimize \ 2 ) with respect to the remaining parameters from the one that 
one want to determine. Since the allowed region in sin 2 2# 13 — p space is already rather 
compact with use of the data of only a far detector the method may produce a satisfactory 
result by just marginalizing the \ 2 . 

Nonetheless, we describe here an alternative iterative method for determination of earth 
matter density in neutrino factory. We start from a zeroth-order assumption of the matter 
density taken from estimate (with uncertainties) by geophysical models. Then, the iterative 
procedure includes the following two steps: 

• We assume an n-th order assumption of the matter density with uncertainties and 
carry out the analysis to obtain n-th order values of #13 and S. 

• We use n-th order values of #13 and 5 with uncertainties to obtain (n + l)-th order 
value of the matter density. 

Assuming convergence we expect that the method is able to produce an improved estimate 
of the three relevant parameters, 8 13 , 5 and the matter density p in a self-consistent way. 

We note that, irrespective of which analysis method is employed, one can expect im- 
provement of the accuracies in determination of the matter density. Measurement of the 
intermediate detector should improve Sp/p directly because of its sensitivity to the earth 
matter effect, in a similar manner as CP phase measurement by the two-detector method 
I32{|. 6 We have already mentioned in the previous section the reasons why a significant im- 



provement can be expected at small # 13 by combined analysis where a intricate structure of 
X 2 minima entails worsened errors in a limited region of 5. 



See Fig. 2 in [53| obtained in the analysis for the first reference in 39( for demonstration of this point 
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IX. CONCLUDING REMARKS 



In this paper, we have demonstrated that an in situ precision determination of the earth 
matter density can be carried out by neutrino factory with a detector placed at the magic 
baseline. There are two regions which are separated by a critical value of #13, sin 2 2# 13 ~ 
a few x 10" 3 . In the large (9 13 region, the uncertainty 5p/p has only mild (^-dependence as 
expected by the desirable feature of measurement at the magic baseline. The achievable 
accuracy of matter density determination is excellent; The fractional error Sp/p is about 
0.43%, 1.3%, and < 3% at la CL at sin 2 29 13 = 0.1, 10" 2 , and 3 x 10~ 3 , respectively. It 
is worth to note that the sensitivity is insensitive against increasing the systematic errors 
by a factor of 2. These are the good enough accuracies to solve the problem of notorious 
uncertainty in the neutrino factory data analysis for precision determination of #13 and the 



CP phase S in view of Fig. 24 of [27[ , for example. 

In the smaller # 13 region, however, we uncovered a subtle (^-dependence in the accuracy 
of determination at small #13. In a certain range of 5 it occurs that the responses of the 
dominant atmospheric term and the solar-atmospheric interference term to matter density 
variation cancel with each other, leading to reduced sensitivity (spike structure in Fig. H]) 
in Sp/p. Adding the disappearance channel does not appear to help. We were unable 
to solve this problem within the framework of analysis used in this paper. We, however, 
emphasize that 5p/p remains small in most of the region of <5; At sin 2 2#i 3 = 10 -4 it is 
3%-7% at la CL in more than 3/4 of the entire region of 5. Furthermore, we expect that 
simultaneous analysis of the intermediate and the far detector data should improve the 
situation significantly, as discussed in Sec. IVHI If it works it resolves the long debated 
problem of obscured CP violation by the uncertainty of the earth matter density. 

It should be stressed, among other things, is that this apparatus seems to provide the 
most accurate way to directly measure the matter density inside the earth. It will provide 
a stringent test for geophysical estimation of the matter density and/or the MSW theory of 
neutrino propagation in matter. 

Finally, we want to mention about other possible ways which could provide with us 
different ways of measuring matter density or alternative means of testing the MSW theory. 

• The transition region of 8 B neutrinos from the sun in the energy region of 3 — 4MeV, if 
measured with precision, would allow us to accurately determine the matter density in 
the solar interior. Given the fact that we have two alternative ways to cross check the 
results, helioseismology and the standard solar model calculation, this possibility may 
give another good way for stringent test of the MSW theory of neutrino propagation 
in matter. 

• As one can learn from the expression of the oscillation probability ffT2l) and ({TBI) the 
relative importance of the matter effect to the vacuum effect in neutrino oscillation 
is not controlled completely by the ratio aL/A 3 i. It may open the possibility of 
using other variables as tunable parameters, e.g., the baseline L. It would then be 
interesting to examine this possibility in the light of recent proposal of the Tokai-to- 



Kamioka-Korea (T2KK) project [39|, [54j. Because of the identical two detector setting, 
it would allow us to accurately determine the matter density in the crust region below 
the sea of Japan. 

Finally, we want to note that measuring the matter density and confronting it to other 
method is not the only way to verify the MSW theory. For example, the concept of 
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mass eigenstate in matter as well as the regeneration of neutrino flavors are best tested 
by observing the day-night effect in solar neutrino observation. 
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